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Abstract 

We present a study of several systems in which a large scale field is generated 
over a turbulent background. These large scale fields usually break a symmetry 
of the forcing by selecting a direction. Under certain conditions, the large scale 
field displays reversals so that the symmetry of the forcing is recovered statistically. 
We present examples of such dynamics in the context of the dynamo instability, 
of two dimensional turbulent Kolmogorov flows and of turbulent Rayleigh-Benard 
convection. In these systems reversals occur respectively for the dynamo magnetic 
field, for the large scale circulation generated by a periodic forcing in space and for 
the large scale roll generated by turbulent thermal convection. We compare the 
mechanisms involved and show that their properties depend on some symmetries 
of the system and on the way they are broken. 
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I. INTRODUCTION 



It is often believed that fully developed turbulent flows are statistically invariant under 
symmetry transformations of the forcing that generates them. Indeed, although the tran- 
sition to turbulence involves successive bifurcations, each of them breaking some spatial or 
temporal symmetry, it is often observed that these symmetries are recovered in a statistical 
sense when the Reynolds number is increased further. The common belief is that strong 
turbulent fluctuations trigger random transitions between symmetric solutions and thus pre- 
vent the flow from staying in a regime with a broken symmetry. In other words, a strongly 
turbulent flow is expected to explore all the available phase space. However, it has been 
known since a long time that this is not so simple and that clear-cut transitions keep occur- 
ring within the strong turbulent regime. The oldest example is provided by the drag crisis 
(see for instance Tritton 1977). The mean drag of a sphere or a cylinder in a turbulent flow 
drops suddenly for a critical value of the Reynolds number Re of order 10^ (the pre-factor 
depending on surface properties) . This corresponds to a transition where the mean flow pat- 
tern changes, the wake becoming narrower. Another example is related to Rayleigh-Benard 
convection, i.e. the flow generated by heating from below a horizontal layer of fluid. It has 
been observed that for a Rayleigh number of order 10^, i.e. roughly 1000 times larger than 
the critical Rayleigh number for the onset of convection, a large scale flow is generated with 
a horizontal extension equal to the length of the container (Krishnamurti & Howard 1981). 
Von Karman swirling flows, i.e. flows generated in a cylindrical volume by the rotation of 
two co-axial disks, also display turbulent bifurcations. In the case of co-rotating disks, an 
axisymmetric mean flow with a strong axial vortex is observed. When the rotation rates are 
varied, this flow breaks axisymmetry, thus generating a roughly periodic modulation of the 
turbulent velocity field superimposed to turbulent fiuctuations with Re ~ 10^ (Labbe et al. 
1996). In the case of disks counter- rotating at the same frequency, the forcing is symmetric 
with respect to a rotation of angle tt about any radial axis in the mid-plane between the two 
disks (see below) . It has been found that this symmetry can be broken through a bifurcation 
that occurs for Reynolds number in the range 10^ to 10^ (Ravelet et al. 2004). 

There are a lot of other examples of transitions leading to broken symmetries in strongly 
turbulent fiows. Bifurcations from a turbulent regime are therefore commonly observed. 
However, in contrast to bifurcations from stationary or space and time periodic flows, that 
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are well documented and for which well-known techniques exist to handle both the linear 
stability problem and the weakly nonlinear bifurcated regime, the concept of bifurcation from 
a fully developed turbulent flow is more questionable even at the level of a clear cut definition. 
In this respect, the dynamo effect in a liquid metal provides a very interesting example in 
which a magnetic field is generated by an instability process that occurs when the kinetic 
Reynolds number of the flow is usually larger than 10^. Although the experiments involve 
some cost and technical difficulties, once the dynamo regime is reached, the dynamics of the 
magnetic field can be easily measured. Then, several interesting problems can be addressed. 
One of them is related to the possibility of anomalous scalings of the magnetic energy density 
above the dynamo threshold due to turbulent fluctuations (Petrelis et al. 2007). Another 
one concerns the effect of turbulent fluctuations on the dynamics of large scale magnetic 
field. The symmetry B(r,t) — B(r, of the equations of magnetohydrodynamics is 
spontaneously broken at the dynamo onset. Slightly above the threshold of a supercritically 
bifurcating dynamo, the magnetic energy density is much smaller than the kinetic energy 
density of turbulent fluctuations. However, turbulent fluctuations are not able to trigger 
a direct transition from B to —B as would be observed in the naive analogy of a particle 
in a two-well potential in the presence of noise. In contrast, the VKS experiment has 
shown that the reversal trajectories display very robust features that are not smeared out 
by turbulent fluctuations (Berhanu et al. 2007). In other words, the dynamics of the large 
scale magnetic fleld involves only a few modes that look weakly coupled to the turbulent 
background, i.e. the reversal dynamics takes place in a low dimensional phase space. Thus, 
the VKS dynamo provides an example in which a few magnetic modes are governed by a 
low dimensional dynamical system although they occur on a strongly turbulent background 
at kinetic Reynolds number larger than 10^. 

Similar features are also observed on purely hydrodynamical models and the same prob- 
lems can be addressed. This is the case in many geophysical and astrophysical flows where 
an oscillatory or a weakly chaotic behaviour can occur in flows at huge Reynolds numbers. 
Some climatic phenomena indeed display a characteristic feature of low dimensional chaos: 
well defined patterns occur within a random temporal behaviour. Examples arc atmospheric 
blockings that can affect the climate on a time scale of several days (Ghil & Childress 1987) 
or El Nino events that occur every few years (VaUis 1986). The qualitative features of these 
phenomena have been often modelled using a few coupled variables such as mean temper- 
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ature, wind or current. The quasi-biennial oscillation(QBO) provides a striking example 
of a large scale almost cyclic reversing flow in the otherwise turbulent atmosphere. It is 
a roughly periodic oscillation in the strength and direction of the zonal (east-west) wind 
in the lower and middle stratosphere over the equator of the Earth's atmosphere. It has 
been observed for more than 50 years in climatological records. Its period fluctuates, the 
mean being slightly larger than 2 years. A reversal flrst appears at an altitude of roughy 
40 km and then propagates downward at a rate of 1 km/month. Thus, these reversals 
are related to a downward drifting pattern. The QBO has been understood to arise from 
the interaction between upward-propagating waves, generated in the troposphere, and the 
mean zonal flow at upper levels where the waves are dissipated (for a review, see Baldwin 
et al. 2001). Reversals of large scale flows driven on a turbulent background, have been 
also observed in laboratory experiments, such as thermal convection (see Ahlers et al. 2009 
for a review) and Kolmogorov flows, i.e., quasi-two-dimensional flows generated by forcing 
an array of counter-rotating vortices (Sommeria 1986). We flrst recall here old and recent 
observations related to bifurcations generating reversals of a magnetic fleld or of a mean flow 
in the presence of turbulent fluctuations and we discuss them in the framework presented 
above. Experimental results and models for reversals of the magnetic fleld are presented in 
section [Tll Reversals of the large scale flow generated by Kolmogorov forcing are studied in 
section IIIII An experimental conflrmation of Sommeria's results is provided and compared 
with a direct numerical simulation. A description of these reversals using a low dimensional 
dynamical system is given and their main features are compared to the ones observed for 
the magnetic fleld. Reversals of the large scale flow in Rayleigh-Benard convection are con- 
sidered in section IIVI After recalling the results obtained in the literature, we present a new 
experiment showing that three-dimensional perturbations can play an important role in the 
reversal process. Finally, we discuss the problem of reversals in relation to the phenomenon 
of drifting patterns and using the framework of amplitude equations. Thus, different types 
of bifurcations leading to reversals in this variety of systems can be easily identifled. We 
conclude with a discussion on the relative contributions of the deterministic dynamics and 
of the turbulent fluctuations to reversals observed in strongly turbulent flows. 
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II. REVERSALS OF THE MAGNETIC FIELD GENERATED IN A FLUID DY- 
NAMO 



A. Dynamos and reversals of the magnetic field 

The dynamo instability is a generic mechanism that converts motion, i.e. kinetic energy 
into electric current, i.e. magnetic energy. This process has been used for a long time to 
generate electricity from the motion of solid rotators. Larmor (1919) proposed that in the 
Sun a similar instability could be operating that would convert kinetic energy of the flow of 
an electrically conducting fluid into magnetic energy. Nowadays it is strongly believed that 
most astrophysical objects (planets, stars, galaxies, ...) generate a magnetic fleld through a 
dynamo instability (Moffatt 1978). In the last decade, the flrst three experimental observa- 
tions of fluid dynamos have been achieved in turbulent flows of liquid sodium (Gailitis et al. 
2001, Stieglitz and Muller 2001, Monchaux et al. 2007). 

Since the work of Brunhes (1906), it is known that the Earth magnetic fleld keeps a 
roughly constant direction for long durations but from time to time it reverses and evolves 
toward the opposite direction. Reversals of the Earth magnetic fleld have motivated a lot 
of studies ranging from paleomagnetism to numerical simulations of models of the Earth 
interior (see for instance Dormy et al. 2000, Petrelis and Fauve 2010). 

Reversals of the magnetic fleld generated by a fluid dynamo have been observed only 
recently in laboratory experiments (the VKS experiment, Berhanu et al. 2007). This ex- 
periment involves a turbulent swirling flow of liquid sodium, generated by two soft iron 
impellers, counter-rotating at frequency Fi (respectively F2) in an inner copper cylinder, as 
sketched in flg. [TJ When the disks counter-rotate with the same frequency F, a statistically 
stationary magnetic fleld is generated when F is large enough. Its mean value involves a 
dominant poloidal dipolar component. Bp, along the axis of rotation, together with a re- 
lated azimuthal component Be, as displayed in flg. [T](left). When the rotation frequencies 
are different, the magnetic fleld can display periodic or random reversals. A time series 
of the fleld in the regime of random reversals is displayed in flg. [2l The fleld fluctuates 
around some value for long durations, typically hundreds of seconds, and suddenly reverses 
within a few seconds. Although the kinetic Reynolds number is large. Re ~ 5 10® for these 
flows, it has been observed that the low dimensional dynamics of the magnetic fleld is not 
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FIG. 1 Sketch of the large scales of the eigenmodes of the VKS experiment. The two disks counter- 
rotate with frequency -Fi and F2. Left: dipolar part of the magnetic mode. Right: quadrupolar 
part. Poloidal, Bp^ and toroidal, Bq^ components are sketched. 

smeared out by strong turbulent fluctuations. In particular, all reversals involve the same 
transitional fleld morphology: the amplitude of the dipolar fleld flrst decreases. If it changes 
polarity, the amplitude increases on a faster time scale and then displays an overshoot before 
reaching its statistically stationary state. Otherwise, the magnetic fleld grows again with 
its direction unchanged. The above features are also observed in recordings of the Earth's 
magnetic fleld, the aborted reversals being often called excursions (Valet et al. 2005). 

B. A generic mechanism for reversals of the magnetic field 

The most striking feature of the VKS experiment is that time dependent magnetic flelds 
are generated only when the impellers rotate at different frequencies (Berhanu et al. 2007, 
Ravelet et al. 2008). We have shown in Petrelis & Fauve (2008) that this is related to the 
broken invariance under TZ-^ when Fi ^ F2 (rotation of an angle vr along any radial axis in 
the mid-plane). In that case, symmetric and antisymmetric modes (under T^^r) are linearly 
coupled. Such modes are displayed in flgure [TJ a dipolar mode is changed to its opposite by 
TZt^, whereas a quadrupolar mode is unchanged. 

Although the flow in the VKS experiment strongly differs from the one in the Earth's core, 
dipolar and quadrupolar modes can be deflned in both cases (using different symmetries, see 
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FIG. 2 Left: Time series of the three components of the magnetic field measured in the VKS 
experiment in the regime of reversals. The two propellers counterrotate with different speeds. 
Right: zoom on two reversals and an excursion (data from Berhanu et al. 2007). 

Petrelis et al. 2009). We assume that the magnetic field is the sum of a dipolar component 
with an amplitude D and a quadrupolar one, Q. We define A = D + iQ and we assume that 
an expansion in power of A and its complex conjugate A is pertinent close to threshold in 
order to obtain an evolution equation for both modes. Taking into account the invariance 
B — )■ — B, i.e. A ^ —A, we obtain 

A = fiA + uA + f3iA^ + f32A^A + I3^AA^ + l3iA' , (1) 

where we limit the expansion to the lowest order nonlinearities. In the general case, the 
coefficients are complex and depend on the experimental parameters. 

Symmetry of the experiment with respect to TZ-^ amounts to constraints on the coefficients. 
Applying the transformation TZ-,^ changes D and Q in different ways: D — —D^ Q ^ Q, 
thus A — )■ —A. We conclude that, in the case of exact counter-rotation all the coefficients 
are real. More generally, the real parts are even and the imaginary parts are odd functions 
of the frequency difference f = Fi — F2. 

The coefficients of ([I]) can be chosen so that it has two stable dipolar solutions ±D and 
two unstable quadrupolar solutions ±Q when / = 0. When / is increased, these solutions 
become more and more mixed due to the increase in the strength of the coupling terms 
between dipolar and quadrupolar modes. Dipolar (respectively quadrupolar) solutions get 
a quadrupolar (respectively dipolar) component and give rise to the stable solutions ±Bs 
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FIG. 3 A generic saddle- node bifurcation in a system with the B — > — B invariance: below thresh- 
old, fluctuations can drive the system against its deterministic dynamics (phase a). If the effect of 
fluctuations is large enough, this generates a reversal (phases b and c). Otherwise, an excursion 
occurs (phase a'). 

(respectively unstable solutions ±Bu) displayed in figure [31 When / is increased further, a 
saddle-node bifurcation can occur, i.e. the stable and unstable solutions collide by pairs and 
disappear (Petrelis & Fauve 2008). A limit cycle is generated that connects the collision 
point. Be, with its opposite. 

This provides an elementary mechanism for field reversals in the vicinity of a saddle-node 
bifurcation. First, in the absence of fluctuations, the limit cycle generated at the saddle- 
node bifurcation connects ±Bc. This corresponds to periodic reversals. Slightly above the 
bifurcation threshold, the system spends most of the time close to the two states of opposite 
polarity ±Bc. Second, in the presence of fluctuations, random reversals can be obtained 
slightly below the saddle- node bifurcation. Bu being very close to Bs, even a fluctuation 
of small intensity can drive the system to B^ from which it can be attracted by — -Bs, thus 
generating a reversal. 

The effect of turbulent fluctuations on the dynamics of the two magnetic modes governed 
by (fl]) can be easily modelled by adding some noisy component to the coefficients. Random 
reversals are displayed in figure S] (left). The system spends most of the time close to the 
stable fixed points ±Bs. We observe in figure H] (right) that a reversal consists of two phases. 
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FIG. 4 (left) Reversals of the magnetic field modelled by ([T|) with an extra noise term that takes 

into account the effect of the turbulent fluctuations, (right) Zoom on reversals and excursions. 

In the first phase, the system evolves from the stable point Bg to the unstable point Bu (in 
the phase space sketched in figure [3]). The deterministic part of the dynamics acts against 
this evolution and the fluctuations are the motor of the dynamics. This phase is thus slow. 
In the second phase, the system evolves from B^ to —Bg, the deterministic part of the 
dynamics drives the system and this phase is faster. 

The behaviour of the system close to Bg depends on the local flow in phase space. Close 
to the saddle-node bifurcation, the position of Bg and B^ deflnes the slow direction of the 
dynamics. If a component of Bu is smaller than the corresponding one of Bg, that component 
displays an overshoot at the end of a reversal. In the opposite case, that component will 
increase at the beginning of a reversal. For instance, in the phase space sketched in flgure 
[3l the component D decreases at the end of a reversal and the signal displays an overshoot. 
The component Q increases just before a reversal. 

For some fluctuations, the second phase does not connect Bu to —Bg but to Bg. It is 
an aborted reversal or an excursion in the context of the geodynamo. Note that during the 
initial phase, a reversal and an excursion are identical. In the second phase, the approaches 
to the flxed point differ because the trajectory that links Bu and Bg is different form the 
trajectory that links Bu and —Bg. In the case of flgure [3l the dipole displays an overshoot 
at the end of a reversal and reaches smaller values during an excursion (see flgure H] right). 
In contrast the quadrupole displays an overshoot at the beginning of a reversal and reaches 
larger values during an excursion. 
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The above mechanism explains several observations made in numerical simulations: a 
dipole-quadrupole interaction is clearly visible in the first reversals simulated by Glatzmaier 
& Roberts (1995). They notice that "the toroidal field is asymmetric with respect to the 
Equator before and after the reversal but is symmetric midway through the transition" . The 
field has thus has a quadrupolar symmetry at the transition (see their figure 2). This has 
been confirmed by the simulations of Sarson & Jones (1999) who find that reversals rely 
"heavily upon the interaction between dipole and quadrupole symmetries" and that they 
are triggered by the random emission of poleward light plumes, i.e. events that break the 
equatorial symmetry of the flow. Similar features have been observed by Wicht & Olson 
(2004). Li et al. (2002) also emphasise that "the dipole polarity can reverse only ... where 
the north-south symmetry of the convection pattern is broken" and that "the quadrupole 
mode grows ... before the reversal". It has been also shown that if the flow or the magnetic 
fleld is forced to remain equatorially symmetric, then reversals do not occur (Nishikawa & 
Kusano 2008). A direct numerical simulation of the dynamo generated by a flow driven in a 
sphere by two counter-rotating co-axial propellers (Gissinger et al. 2010) has also displayed 
the following features. Reversals of the axial dipole occur only when the propellers are 
rotated at different rates. When the magnetic Prandtl number {Pm) is small enough, they 
involve an axial quadrupole and the dynamics of the dipolar and quadrupolar modes during 
a reversal is similar to the one observed in the VKS experiment. 

C. A dynamical system displaying random reversals of the magnetic field 

The three modes that are important in the low Pm simulations are the dipole D, the 

quadrupole Q, and the zonal velocity mode V that breaks the TZj^ symmetry. These modes 
transform as D ^ — D, Q ^ Q and V — )■ —V under the TZ^^ symmetry. Keeping nonlinear 
terms up to quadratic order, we get the dynamical system (Gissinger et al. 2010) 

b ^ HD- VQ, 
Q = -uQ + VD, 

V = r-V + QD. (2) 

A non zero value of F is related to a forcing that breaks the TZ^^ symmetry, i.e. propellers 
rotating at different speeds. 
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FIG. 5 Numerical integration of the amplitude equations (2). Time recording of tlie amplitude of 
the quadrupolar mode for /i = 0.119, v = 0.1 and T = 0.9. (Figure from Gissinger et al. 2010). 

The dynamical system (2) with F = occurs in different hydrodynamic problems, such 
as the dynamics of periodic structures in thermoaline convection. It has been analysed in 
detail (Hughes & Proctor 1990). The relative signs of the coefficients of the nonlinear terms 
have been chosen so that the solutions do not diverge when fi > and z/ < 0. Their modulus 
can be taken equal to one by appropriate scalings of the amplitudes. The velocity mode is 
linearly damped and its linear coefficient can be taken equal to —1 by an appropriate choice 
of the time scale. Note that similar equations were obtained through drastic truncation of 
MHD equations (Nozieres 1978). However, in that context fi should be negative and the 
damping of the velocity mode is discarded, thus strongly modifying the dynamics. 

This system displays reversals of the magnetic modes D and Q for a wide range of 
parameters. A time recording is shown in figure [5l The mechanism for these reversals 
results from the interaction between the modes D and Q coupled by the broken 7^^ symmetry 
when V ^ 0. It is thus similar to the one described above but also involves an important 
difference: keeping the damped velocity mode into the system generates chaotic fluctuations. 
It is thus not necessary to add external noise to obtain random reversals. This system is 
fully deterministic as opposed to the one of Petrelis and Fauve (2008). Its phase space shows 
the existence of chaotic attractors in the vicinity of the ±B quasi-stationary states. When 
these symmetric attractors are disjoint, the magnetic field fluctuates in the vicinity of one 
of the two states ±B and the dynamo is statistically stationary. When /x is varied, these 




11 



two attractors can get connected through a crisis mechanism (Grebogi et al. 1982), thus 
generating a regime with random reversals. Compared to previous deterministic models, 
the one obtained here displays dynamical and statistical properties that are much closer 
to the ones of the direct simulations at low of Gissingcr et al. (2010) and of the VKS 
experiment. For instance, the direct recordings oi D oi Q do not involve the growing 
oscillations characteristic of reversals displayed by the Rikitake (1958) or Lorenz (1963) 
systems but absent in dynamo experiments or in direct simulations. Correspondingly, the 
probability density function of D is also much closer to the one obtained in experiments or 
direct simulations than the one of previous deterministic models. 



III. LARGE SCALE CIRCULATION OVER A TWO DIMENSIONAL TURBU- 
LENT FLOW 

A. Large scale quasi-two-dimensional flows 

Rapid rotation or an externally applied strong magnetic field have been known since a long 
time to make the flow two-dimensional by inhibiting velocity gradients along the direction 
of the rotation axis or of the magnetic held (Chandrasekhar 1961). This situation is of 
obvious interest for geophysical or astrophysical flows for which rotation can be important. 
The properties of turbulent flows strongly differ between these nearly two dimensional (2D) 
systems and the usual three dimensional ones. In 2D, it is well known that an inverse 
cascade tends to accumulate energy in the largest scales of the system. Early studies of 
periodically driven flows (Sommeria 1986, Tabehng et al. 1987) have shown that an array 
of vortices becomes unstable when the forcing is increased and that a large scale circulation 
is generated in the whole system. This large scale circulation is expected to display rich 
dynamics because it is coupled to the smaller scale turbulent fluctuations. Both the large 
scale field and the fluctuations are therefore components of the same turbulent fleld, the 
velocity fleld. This differs from the the former example in which reversals involve a fleld 
(the magnetic fleld) that is different from the turbulent velocity fleld. 
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B. Reversals of the large scale circulation in Kolmogorov flows: experiments and 
numerical simulation 

Following Sommeria (1986) and with the purpose of studying the dynamics of the large 
scale circulation, we consider a square layer of mercury of length L = 12 cm and depth 
a = 1.6 cm. This layer is placed into a vertical magnetic field up to i? = 500 G. An electric 
current / up to 40 A is injected in the fluid through 36 electrodes placed on a periodic square 
array. This current generates a Laplace force which drives the motion of the fluid. When the 
applied magnetic field is large enough, the flow is mostly two dimensional. At low current, 
a flow made of 36 counterrotating vortices is generated. This velocity field, obtained from 
direct numerical simulations presented below, is shown in fig. [7] (a). At larger current, a 
large scale flow is generated. Two states are possible and are shown in fig. [7] (b,c). Let 
(x, y) be the coordinates measured from the lower left corner of the square. We note that the 
forcing is invariant under the reflections with respect to the planes x = L/2 and y = L/2. 
We call Sx (respectively Sy) these symmetries. The large scale circulation thus breaks the 
symmetry of the forcing by selecting a direction of rotation. 

For large Reynolds number Re = {BIL^ /(paiy'^)y^'^, there exists a critical value of Rh = 
(al /(Bai'L'^)Y^'^ under which the flow displays apparently random reversals between the 
two possible large scale circulations (Sommeria 1986). We present measurements of the x 
component of the velocity at position (6; 2) cm and (6; 10) cm. The velocity is measured 
using Vives probes and is averaged over 2 cm. In fig. [6l a time series is displayed in the 
regime of random reversals. The direction of the large scale circulation remains constant 
for long durations, here larger than one hour, and reverses in less than 100s. Note that the 
shape of the reversals does not display a clear overshoot. In the case of the VKS experiment, 
reversals are observed when the two disks rotate at different speeds. Then the TZ-,^ symmetry 
is broken and this constrains the form of the phase space (fig. [3]) which enforces the shape of 
the reversals and the existence of an overshoot. In contrast, in the case of the 2D turbulent 
flow, no externally broken symmetry imposes an overshoot. 

A numerical simulation of this experiment has been performed. When the magnetic field 
is strong enough, the flow is almost 2D and the stream function ijj obeys the two-dimensional 
Navier-Stokes equation 

dtAip + J{Aip, ip) = -^Aijj + + 6vr sin(67rx) sin(67ry) . (3) 

Rh Re 
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FIG. 6 Time series of the velocity V jV^ where K = {Bl/ipa))^/"^ ~ 4.28 cm.s"^ for Rh = 37.7 
and Re = 45000. The two curves correspond to measurements located on opposite sides of the cell. 
The velocity is averaged between one centimeter and three centimeters away from the lateral wall. 

The first term on the right hand side is linear friction coming from the Hartmann layers on 
the bottom of the cell. The second one is the viscous term, and the last one mimics the 
electromagnetic forcing generated by the 6x6 electrodes. We solved this equation using a 
pseudo-spectral code with uniform time-stepping and stress free boundary conditions (that is 
vanishing normal velocity and tangential constraint at the boundaries). The stream function 
is decomposed in the domain {x,y) G [0;!]^ on the basis (sin(n7rx) sin(p7ry))j.j^p-jg|-,^ ^p. 
The linear terms are computed in Fourier space, while the nonlinear one is computed in real 
space. All the numerical runs described further were performed using 256 Fourier modes in 
each direction. We checked convergence: doubling this number of modes does not modify 
the numerical solution. 

C. A dynamical system displaying reversals of the large scale circulation 

Reversals of the large scale circulation are observed in this direct numerical simulation. 
An example of time series is shown in fig. |8l It is possible to identify the large scale structures 
that are involved in the process. They can be decomposed into three main components : 
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FIG. 7 Stream function of the numerically computed solution of equation [3] for Re = 2000, Rh = 55, 
a. Snapshot at short time, the velocity field displays the same periodicity as the forcing, b and 
c are two snapshots displaying large scale circulations of opposite signs. The large scale flow has 
broken symmetries of the forcing. Under certain conditions, the circulation reverses from one state 
to the other. 
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FIG. 8 Left: Time series of the amplitude of the large scale circulation D and of two large scale 
components Qx and Qy (see definition in the text) obtained numerically for Rh = 55 and Re = 2000. 
Right: zoom on one reversal. A sliding average is applied to reduce high frequency fluctuations. 

• The amplitude of the large-scale circulation is denoted as D and is defined here as the 
Fourier coefficient of sin(7rx) sin(7r?/) in the decomposition of the streamfunction. 

• A second component, called Q^, is the part of the large scale flow that is quadrupolar 
in the x direction. It is even under Sx and odd under Sy. It can be measured by the 
Fourier coefficient of ip in sin(27ra;) sin(7r?/). 

• A third component, called Qy, is the part of the large scale flow that is quadrupolar 
in the y direction. It is odd under Sx and even under Sy. It can be measured by the 
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Fourier coefficient of in sin(7rx) sin(27ry). 

We observe that reversals of D take place when Qy reaches large values. In agreement with 
the experiment, no overshoot is seen at the end of the reversals. 

In order to obtain the minimal model that displays similar dynamics, we note that the 
three components change differently under the different symmetries of the system. As al- 
ready mentionned D and Qx are odd under Sy while Qy is even. D and Qy are odd under 
Sx while Qx is even. A simple dynamical system that verifies these symmetries is 

ID = -vD- QxQy 

Qx Qx QyD Qx 

Qy = fiQy + DQx. (4) 

Compared to Hughes and Proctor (1990), the set of equations (jl]) contains an extra third- 
order nonlinearity in the last equation. A time series is represented in figure [9] for fi = 6 
and u = 2.65: the large scale circulation D reverses randomly. As in the direct numerical 
simulations, the value of Qy increases during the reversals of D. 

One should notice that the present problem maps onto the magnetic one in the following 
way: if we replace D, Qy, and Qx respectively by the amplitude of the dipolar magnetic 
component D, quadrupolar magnetic component Q, and equatorially-antisymmetric veloc- 
ity component V, the symmetries and Sy lead to the same constraints on the dynamical 
system than the invariance B — t- — B and the equatorial symmetry for the magnetic field of 
an astrophysical object. Therefore the system of equations (jlj) could describe the magnetic 
field of an astrophysical object as well, and leads to reversals of the magnetic dipole with 
no imposed equatorial-symmetry breaking. In this sense it is slightly different from model 
(2) in which the dynamics of V (corresponding to Q^) results from an externally imposed 
broken symmetry and not from an instability of the large scale velocity field that sponta- 
neously breaks this symmetry. In these low dimensional models, the system follows chaotic 
trajectories that wander in the basin of attraction associated to one polarity of the large 
scale field or its opposite. In the regime of reversals, the chaotic attractor connects the two 
basins of attraction and after a long duration the system escapes from one state toward its 
opposite. 
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FIG. 9 Time series of D solution of equations (j3]) in a regime of reversals (/z = 6 and v = 2.65). 

IV. LARGE SCALE CIRCULATION IN TURBULENT RAYLEIGH-BENARD 
CONVECTION 

A. Experiments on the dynamics of the large scale circulation in Rayleigh-Benard 
convection 

Rayleigh-Benard convection is a canonical example of cellular instability. It is achieved 
by uniformly heating from below a horizontal layer of fluid. For small temperature differ- 
ence AT across the layer of depth d, the fluid remains in a stable heat conducting state, 
with a linear temperature profile and no fluid motion. However, if the fluid's density is a 
decreasing function of temperature, the thermal gradient generates a density stratification 
with cold heavy fluid above warm light fluid. For sufficiently large temperature differences, 
the resulting buoyancy force overcomes dissipative effects due to kinematic viscosity, z/, and 
heat diffusivity, k, causing less dense warmer fluid to rise and cooler fluid to sink. When the 
top and bottom boundaries have a thermal conductivity much higher than the one of the 
fluid, i.e. their temperatures are homogenous, periodic parallel convection rolls of horizontal 
size comparable to rf, result from the circulation of the fluid, when the Rayleigh number, 
Ra = agATd^/luK,), is larger than a critical value, Rac {g is the acceleration of gravity, a is 
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the coefficient of thermal expansion of the fluid). In the case of a large aspect ratio container 
(ratio of the horizontal size L to the depth ci), it has been observed by Krishnamurti and 
Howard (1981) that a large scale horizontal shear flow of horizontal size L is generated for 
large enough values of Ra, roughly three orders of magnitude larger than Ra^.. In contrast 
to the large scale flow generated by the Kolmogorov forcing considered in section llllt the 
dynamics of this flow has not been studied. For convenience, most high Rayleigh number 
experiments have been performed in cylindrical containers of aspect ratio one. In that case, 
the large scale circulation evolves out of the cellular structure observed at convection onset 
without implying a secondary bifurcation and an associated broken symmetry. However, 
even in that simple case, this large scale circulation displays dynamics that have been first 
related to flow reversals (see Ahlers et al. 2009 for a review of the early experiments and 
related models). It has been realised that, although simpler for the machine shop, the cylin- 
drical geometry is not appropriate to study flow reversals. The convective pattern at onset 
and thus the large scale circulation for large Ra indeed break the rotational invariance of the 
container. Thus, the dynamics primarily results from the azimuthal orientation of the large 
scale circulation and the flow measured locally can vanish without necessarily implying a 
cessation of the whole large scale circulation (Ahlers et al. 2009). Both reorientations and 
cessations of the large scale circulation have been evidenced by recent direct numerical sim- 
ulations (Mishra et al. 2010). They also displayed a quadrupolar temperature distribution 
during flow cessations (as opposed to the dipolar distribution resulting from the large scale 
circulation). Thus, although different mechanisms are superimposed in the dynamics of large 
scale circulation in cylindrical containers, flow cessations involve the competition between 
two large scale modes with different symmetries, in striking analogy with the mechanism of 
reversals of the magnetic field proposed in III.BI 

Only a few studies of the dynamics of the large scale circulation have been performed in 
parallelepipedic containers. Liu and Zhang (2008) reported a large scale flow in an elongated 
box of aspect ratio 2.6 for Ra > 10^ together with reversals but they provided a detailed 
study only with freely moving nylon spheres within the flow. Another experiment has 
been performed recently in a quasi-two-dimensional parallelepipedic container of aspect ratio 
one, and its results have been compared to direct numerical simulations of two-dimensional 
convection (Sugiyama et al. 2010). Although the authors insist on the importance of corner 
flows superimposed to the large scale diagonal circulation, a more pertinent observation 
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is related to the change of symmetry of the large scale flow observed during the reversal 
process. When the dipolar part of the large scale flow vanishes, a quadrupolar structure 
is clearly visible (see their figure lb). We also note that earher two-dimensional numerical 
simulations (Breuer & Hansen 2009) performed in a container of aspect ratio two, also 
displayed the competition between two large scale modes with different symmetries during 
the reversal process. Thus, it appears that, as for reversals of the magnetic field in the 
VKS experiment, reversals of the large scale convective flow in the quasi- two-dimensional 
configuration also involve the competition between two large scale modes with different 
symmetries. Another similarity between the convective and dynamo reversals is related to 
their occurrence in a finite interval range of the control parameters (Rayleigh and Prandtl 
numbers Ra, Pr or magnetic Reynolds numbers Rm)- In both two-dimensional numerical 
simulations of convection, it has been also reported that relatively small changes in the 
aspect ratio strongly modify the reversal dynamics and thresholds. This probably results 
from the inhibition of the mode with quadrupolar symmetry that is involved in the reversal 
process. Finally, the numerical simulations of Sugiyama et al. (2010) show that reversals 
disappear when the Prandtl number of the fluid is too small (Pr < 0.7) or too large. We will 
report below reversals in a cubic container of mercury (Pr ~ 0.025). Applying a horizontal 
magnetic field enables to make the flow more and more two-dimensional by inhibiting velocity 
gradients along the direction of the magnetic field (Fauve et al. 1984). We will show that 
this strongly affects the modes involved in the reversal process. 

B. Dynamics of the large scale circulation in a cubic container subject to a hori- 
zontal magnetic field 

We consider a cubic container of size d = 45 mm filled with mercury of den- 
sity p — 13.510^ kg. m~^, kinetic viscosity u — 1.12 10~^m^.s~^, thermal diffusivity 
K, — 4.3 10~^m^.s~^, specific heat Cp — 13SJ.kg~^.K~^ and coefficient of thermal expansion 
a — l.S10~'^K~^. The side walls of the experiment are made of PVC while the top 
and bottom ones are copper plates. Water circulation through the upper copper plate 
maintains the upper temperature fixed. A constant heat flux F is imposed at the bottom 
boundary. The lower plate is covered with temperature probes which allows to determine 
the temperature difference between the bottom and the top plate AT and also to identify 
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the spatial structure of the temperature field in the experiment. A horizontal magnetic 
field, B, up to 500 G can be applied perpendicular to one of the sides of the cell. In fig. [TOl 
we display the Nusselt number Nu = F/ {pCpnATd) as a function of the Rayleigh number. 
In this regime of parameters, the magnetic field has little effect on the heat transport. For 
AT ~ 48°C i.e. Ra = 1.58 10'', the variation of the Nusselt number is less than 4% when B 
is increased from to 500 G corresponding to an interaction parameter, ~ \J ag^^T i 
smaller than 0.2. 
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FIG. 10 a. Nusselt number as a function of the Rayleigh number for (o): B = OG, (★): B = 100 G, 
(□): B = 200 G, (*): B varies from to 500 G at AT ~ 48"G, (o): B varies from to 325 G at 
AT ~ 35°G. b. Parameter space {B,AT). The blue squares correspond to stationary regimes. 
The red circles correspond to reversals of the large scale circulation. 

In the absence of magnetic field, for Ra smaller than Rar — 4.3 10^, the temperatures 
fiuctuate around some constant value. In this regime, the larger scale of the convective 
pattern is a single roll that fills the whole cell. This roll is expected to be perpendicular to 
one of the sides of the cell and we have checked that its direction can be rotated by tilting 
the cell. However, no spontaneous change in the large scales of the temperature pattern is 
observed even for measurement durations larger than 24 hours. 

A change of behaviour is observed when Ra reaches Rar. The temperatures fiuctuate 
for long duration around some value but from time to time, the whole temperature pattern 
changes suddenly. Two time series of the temperature measured at opposite positions on the 
diagonal of the cell, 40 mm away from the corners, are displayed in fig. [TTJ These large scale 
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modifications are separated by very long time intervals, so that statistical properties can only 
be estimated. We have measured that the mean duration between these events decreases 
from more than 3 hrs to around 30 minutes when AT changes from 13.4 to 27.7° C. 




t(hr) 

FIG. 11 Time series of the temperature measured at two opposite corners of the lower plate for 
B = 0, AT = 15.4°C and Ra = 5.1310^. Measurement rate is 1/3 Hz and a sliding average over 
10 measurements is performed to reduce fluctuations. 

Although the time series share common features with the ones of the reversals observed 
in the VKS experiment (see fig. [2]) or in the two dimensional turbulent Kolmogorov flow (see 
fig. [6]), we note that the evolution is more complicated: there are more than two stationary 
states and we can identify events where one of the temperatures changes with no noticeable 
evolution of the other one. Careful inspection of the evolution of several probes shows that 
the modifications of the large scale of the flow are associated to a rotation of the large scale 
roll. This behaviour can be identified by considering the horizontal temperature gradients. 
In fig. [121 we plot a cut in phase space (T^ — T^, T^ — T^) where T^ — Tb and T^ — T^ are 
temperature differences measured along two perpendicular sides of the cell. The trajectories 
in phase space are roughly located along a circle. Between the changes of direction, the 
system fluctuates close to one of the four states that correspond to the roll rotation vector 
pointing towards one of the four lateral sides of the cell. Indeed if a roll is aligned with the 
line a — b, then T^ — T^ is roughly zero and the direction of rotation of the roll is given by 
the sign of T^ — T^. The trajectory in phase space shows that the roll evolves from being 
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perpendicular to one of the sides to being perpendicular to one of the two neighbouring sides. 
Its direction thus changes by ±90°. In other words, the large scale flow displays random 
changes of direction and reversals are obtained by the sequel of two changes of orientation 
in the same direction, which results in a rotation of 180°. 
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b. Ta-Tb (C) 

FIG. 12 Projection of the phase space obtained from temperature differences: Ta — and — T^. 
a- In the absence of applied magnetic field, same data as fig. [TT] b- For AT = 48°C, a magnetic 
field B = 375 G is applied parallel to the direction a — b. 

The obtained phase space clearly displays the four- fold symmetry of the experiment. A 
way to lower the degree of symmetry is to apply a horizontal magnetic field. Measurements 
have been performed in which a field is applied parallel to the line joining a and b. In 
the parameter space displayed in fig. [10] b, we observe that reversals disappear when the 
magnetic field is large enough. When reversals are observed in the presence of an applied 
magnetic field, the trajectories in phase space are modified. As displayed in fig. [12] b, the 
two states corresponding to rolls perpendicular to the field (T^ — T^ — 0) have lost stability: 
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b. 

FIG. 13 Phase space displaying stable solutions (□ black) and unstable solutions (o red), a- The 
four-fold symmetry is the one of the experiment in the absence of magnetic field, b- In the case of 
an applied magnetic field, the two formerly stable solutions representing rolls perpendicular to the 
field are unstable and two pairs of new stable solutions (the diagonal rolls) are formed through a 
pitchfork bifurcation. Considering that turbulent fluctuations allow to evolve from a stable solution 
toward its neighbouring ones, both phase spaces represented in fig. [12] can be reconstructed. 

the system spends only a short time in these states. In contrast the system spends long 
durations in the vicinity of four new states. These states correspond to Ta — Tb ~ ±(T(i — T^) 
which in physical space means that a roll is aligned along the diagonal of the square. These 
observations result from the applied magnetic field that reduces the variations of the velocity 
along the direction of the field. The large scale dynamics is thus modified by the following two 
effects: first, the rolls are preferably aligned along the magnetic field and rolls perpendicular 
to the field are expected to loose their stability; second, turbulence is reduced because 
fluctuations along the direction of the field are damped. These two effects explain that the 
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pair of solutions perpendicular to the field turns unstable and ultimately that large scale 
dynamics is suppressed. 

From the point of view of low dimensional modelling, the regime with B = corresponds 
to a system with four stable fixed points. In a planar model describing the competition 
between perpendicular rolls, these stable points are located between four unstable states 
which are saddles (see fig. [13] a). When a magnetic field is applied, a pair of stable solutions 
(transverse rolls) looses its stability. The four new states (diagonal rolls) that are observed 
can be described as two pairs of stable states that have been generated from a pitchfork 
bifurcation of the transverse rolls. The resulting phase space contains six nodes and, in- 
between, six saddles (see fig. [T3]b). A way to model the random dynamics of the large scale 
fiow is to consider that turbulent fiuctuations act as noise that triggers transitions from one 
state to one of its neighbours. 

As in the case of reversals of the magnetic field, the large scale fiow does not change sign 
directly under the infiuence of turbulent fiuctuations (plumes, etc). This transition always 
involves another (or several other) large scale modes. In the case of square symmetry, 
perpendicular rolls are involved. When the square symmetry is broken by an external 
magnetic field, more complex mixed states are involved in the reversal process. For the 
Prandtl number of mercury, reversals are not observed above a critical value of the magnetic 
field. In contrast to the reversals of the magnetic field in the VKS experiment, an externally 
imposed broken symmetry is not needed to generate reversals of the large scale circulation 
in convection as well as in Kolmogorov fiows. Correspondingly, the reversals do not display 
an overshoot. 

V. DISCUSSION AND CONCLUSION 

We have presented three examples of turbulent fiows in which a field generated at large 
scale reverses. Some properties are common to the three systems. 

In all cases the large scale field breaks a symmetry of the forcing and reversals allow the 
system to statistically recover this symmetry. 

A clear separation of time scales exists between the duration of a reversal and the waiting 
time between reversals. Indeed, random reversals are rare events that are initiated when 
the system wanders sufficiently far away from the basin of attraction of a fixed point. Once 
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this phase has been completed, the evolution is fast because it does not require to reach a 
rarely visited domain in phase space. 

During a reversal the whole field does not vanish. In particular, it is possible to identify 
large scale structures that play a role in the dynamics and are visited by the system. 

The three systems differ from the point of view of symmetries and in the way they are 
broken. The magnetic field and the large scale circulation have two stable states while the 
large scale roll in our convection experiment has four or six stable states. In the case of the 
experimental dynamo, one symmetry of the forcing is externally broken by rotating the two 
disks at different speeds. In the other two systems, one symmetry is spontaneously broken 
and this selects the direction of rotation of the roll's axis or the sign of the transitional 
two cell flow. These symmetry properties are responsible for the differences in shape of the 
reversals observed in the three systems. 

The geometry of the large scale dynamics being governed by a low-dimensional dynamical 
system, we expect that reversals can also occur in laminar flows. This is indeed the case. 
Direct simulations of axisymmetric Rayleigh-Benard convection in a cylindrical geometry 
(Tuckerman & Barkley 1988) as well as the weakly nonlinear analysis of the same system 
(Siggers 2003), have shown that the steady state observed above the convection onset can 
bifurcate to a travelling wave in which the number of rolls oscillates between two adjacent 
values. It has been emphasized that this transition to periodic behaviour is a symmetry- 
restoring bifurcation, as the states in the second half of the hmit cycle are related to those of 
the flrst half by changing the sign of the velocity held (and accordingly of the temperature 
perturbation). A travelling spatially periodic pattern is indeed the simplest way to generate 
reversals of the flow measured at some location. This has been observed in direct numerical 
simulations of Rayleigh-Benard convection with horizontal periodic boundary conditions 
(Paul et al. 2010). In that case, randomly occurring lateral shifts of the roll pattern by half 
a wavelength, lead to global flow reversals of the convective motion. These events are also 
observed in a flnite interval range of Rayleigh number (although Ra is much smaller than 
in the presence of lateral boundaries). At flrst sight, one may consider that the presence 
of lateral boundaries also deeply modifies the qualitative bifurcation mechanism leading to 
fiow reversals but it is worth considering this analogy further. Drift instabilities of cellular 
patterns have been indeed observed in various experiments, such as Couette fiow between 
two horizontal coaxial cylinders with a partly filled gap (Mutabazi et al. 1988), film draining 
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(Rabaud et al. 1990), directional crystal growth (Simon et al. 1988), Faraday waves (Douady 
et al. 1989), etc. They have been understood as a secondary instability that spontaneously 
breaks the reflection symmetry of the pattern which then travels in the direction related to 
this asymmetry (Coullet et al. 1989, Fauve et al. 1991, Caroli et al. 1992). If the reflection 
symmetry, say x — t- — x, is externally broken, then one expects that a spatially periodic 
pattern that breaks the translational invariance along the x-axis, generically drifts along 
the X-axis at onset when the domain is infinite. This is easily understood in the framework 
of amplitude equations that amounts to describe a periodic pattern of wave number k in 
the form A{t) expikx + c.c. + ■ ■ ■ . A{t) is a complex amplitude, c.c stands for the complex 
conjugate of the previous term and the dots represent higher order terms in the vicinity of 
the bifurcation threshold. It is well known that the form of the equation for A in series 
of successive powers of A and A (the complex conjugate of A), results from the broken 
symmetries at the instability onset. If the system is invariant under the transformation 
A — i- —A, this equation, up to leading order nonlinear terms, is similar to ([T]) 

A = ijA + uA + /3iA^ + P2A^A + P:^AA^ + P^A^ , 

where all the coefficients are complex numbers. If the system is translationally invariant 
along the x-axis, the amplitude equation should be invariant by rotation in the complex 
plane, A — )■ Aexpi^, where is an arbitrary phase; thus z/ = /3i = /^s = = 0. If the 
system has parity invariance, x — )■ — x, which implies the invariance A — )■ A for the amplitude 
equation, then all the coefficients should be real. Thus, we get the usual amplitude equation 
that describes the stationary bifurcation of a homogeneous periodic pattern, 

A = + P2M\^ A, (5) 

where all the coefficients are real. When the parity invariance is externally broken, these 
coefficients generically have non zero imaginary parts. When the imaginary part of /i, say 
/Xj, is non zero, then the pattern is generated through a Hopf bifurcation and is thus a 
travelling wave that leads to periodic reversals of the flow at a given location. When the 
lateral boundaries are taken into account, translational invariance is externally broken and 
the leading order additional term that should be considered is vA. This new term indeed 
tends to quench the phase of the pattern. The linear stability analysis of the solution A = 
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gives the dispersion relation for the growth rate s 

- 2/i^,s+ l/ip - |z/|2 = 0, (6) 

where /i^ is the real part of fi. We have a stationary bifurcation for |/i| = |z/| if fir < 0, a 
Hopf bifurcation for /ir = Oif|/i|>|z/|>0, and a co dimension-two bifurcation for /i^ = 
and fif = Writing A = Rexpi(j), the stability of finite amplitude stationary patterns 
can be studied in the phase approximation provided the amplitude R is slaved to the phase. 
The imaginary part of the linear part of ([1]) gives 

= /ij — z/r sin 2(f) + Vi cos 20. (7) 

As described in section III.Bt stationary solutions disappear via a saddle-node bifurcation 
when /i^ = and a limit cycle that corresponds to travelling waves i.e., periodic reversals, 
is generated. This occurs only if /Xj increases faster than Vi when the reflection symmetry 
is externally broken. If |z/j| > the solutions remain stationary. In the context of the 
dipole-quadrupole interaction in the dynamo problem, a broken symmetry that induces 
much larger than has been found as a mechanism for hemispherical dynamos (Gallet & 
Petrelis 2009). In the case of axisymmetric convection, the broken reflection symmetry is 
related to the curvature of the rolls. 

The above phase approximation breaks down in the vicinity of the codimension-two point. 
Another type of bifurcation from stationary solutions to a limit cycle takes place in that 
case (Gambaudo 1985, Guckenheimer & Homes 1986). It is a subcritical Hopf bifurcation 
that can be easily discriminated from the previous scenario; the limit cycle appears with 
a finite period whereas the period diverges when it is generated through a saddle-node 
bifurcation. Stationary and travelling solutions coexist in some parameter range, thus this 
second scenario displays bistability. Both bifurcation types, saddle-node and subcritical 
Hopf, have been reported for the axisymmetric convection problem (Siggers 2003) as well as 
for reversals of the magnetic field in the VKS experiment (Berhanu et al. 2009). In numerical 
simulations of mean-field dynamo models, the occurence of reversals of the magnetic field 
close to a codimension-two bifurcation point has been emphasized by Stefani & Gerbeth 
(2005). 

As mentioned above, drift bifurcations of cellular patterns also occur in a parity invariant 
system provided the reflection symmetry is spontaneously broken through a secondary bifur- 
cation. This completes the analogy between drifting patterns in laminar flows and reversals 
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of large scale fields in the presence of turbulent fiuctuations that can also occur as a result 
of an external or a spontaneous broken symmetry. 

Turbulent fiuctuations thus do not play a major role in the geometry of reversals, i.e. to 
determine the trajectory that connects two solutions of opposite polarity in phase space. In 
all the cases studied here, this connection involves at least two large scale modes. Turbulent 
fiuctuations of course affect the time dependence of reversals. They trigger random reversals 
below the saddle-node bifurcation to a limit cycle. However, very similar random behaviours 
can be obtained with deterministic models involving three large scale modes, so that it would 
not be so easy to determine the source of randomness in a natural system displaying reversals. 

References 

Ahlers, G., Grossmann, S. and Lohse, Heat transfer and large scale dynamics in turbulent 
Rayleigh-Benard convection. Rev. Mod. Phys., 2009, 81, 503-537. 

Baldwin, M. P., Gray, L. J., Dunkerton, T. J., Hamilton, K., Haynes, P. H., Randel, W. J., 
Holon, J. R., Alexander, M. J., Hironta, L, Horinouchi, T., Jones, D. B. A., Kinnersley, J. 
S., Marquardt, C, Sato, K. and Takahashi, M., The Quasi-Biennial Oscillation, Reviews of 
Geophysics, 2001, 39, 179-229. 

Berhanu, M., Monchaux, R., Fauve, S., Mordant, N., Petrelis, F., Chiffaudel, A., Daviaud, 
F., DubruUe, B., Marie, L., Ravelet, F., Bourgoin, M., Odier, Ph., Pinton, J.-F., Volk, R., 
Magnetic field reversals in an experimental turbulent dynamo, Europhys. Lett., 2007, 77, 
59001. 

Berhanu, M., Monchaux, R., Bourgoin, M., Odier, Ph., Pinton, J.-F., Plihon, N., Volk, R., 
Fauve, S., Mordant, N., Petrehs, F., Aumatre, S., Chiffaudel, A., Daviaud, F., DubruUe, B., 
Ravelet, F., Bistability between a stationary and an oscillatory dynamo in a turbulent fiow 
of hquid sodium, J. Fluid Mech., 2009, 641, 217-226. 

Breuer, M. and Hansen, U., Turbulent convection in the zero Reynolds number limit, Euro- 
phys. Lett, 2009, 86, 24004. 

Brunhes, B., Recherches sur la direction d'aimantation des roches volcaniques. J. de Phys. 
Theor. App., 1906, 5, 705-724. 

Caroli, B., Caroli C. and Fauve, S., On the phenomenology of tilted domains in lamellar 



28 



eutectic growth, J. Physique /, 1992, 2, 281-290. 

Chandrasekhar, S., Hydrodynamic and hydromagnetic stability, 1961, (Clarendon Press, Ox- 
ford) . 

CouUet, P., Goldstein, R. E. and Gunaratne, G. H., Parity-breaking transitions of modulated 
patterns in hydrodynamic systems, Phys. Rev. Lett, 1989, 63, 1954-1957. 

Dormy, E., Valet J. -P. & Courtillot V., Numerical models of the geodynamo and observa- 
tional constraints. Geochem. Geophys. Geosyst, 2000, 1, 2 2000GC000062. 

Douady, S., Fauve, S. and Thual, O., Oscillatory phase modulation of parametrically forced 
surface waves, Europhys. Lett, 1989, 10, 309-315. 

Fauve, S., Laroche, C. and Libchaber, A. and Perrin, B., Chaotic phases and magnetic order 

in a convective fluid, Phys. Rev. Letters, 1984, 52, 1774-1777. 

Fauve, S., Douady, S. and Thual, O., Drift instabilities of cellular patterns, J. Physique II, 
1991, 1, 311-322. 

Gallet, B. and Petrelis, F., From reversing to hemispherical dynamos, Phys. Rev. E, 2009, 
80, 035302. 

Gailitis, A., Lielausis, O., Platacis, E., Dement'ev, S., Cifersons, A., Gerbeth, G., Gundrum, 
T., Stefani, F., Christen M. and Will, G., Magnetic field saturation in the Riga dynamo 
experiment, Phys. Rev. Lett, 2001, 86, 3024-3027. 

Gambaudo, J. M., Perturbation of a Hopf bifurcation by an external time-periodic forcing, 
J. Dijf. Eqs., 1985, 57, 172-199. 

Ghil, M. & Childress, S., Topics in geophysical fluid dynamics: atmospheric dynamics, 
dynamo theory, and climate dynamics. Appl. Math. Sci. vol. 60, 1987 (New York: Springer 
Verlag) . 

Gissingcr C, Dormy E., Fauve S., Morphology of field reversals in turbulent dynamos, 
Europhys. Lett, 2010, 90, 49001. 

Glatzmaier, G.A. and Roberts, P.H., A three-dimensional self-consistent computer simula- 
tion of a geomagnetic field reversal. Nature, 1995, 377, 203-209. 

Grebogi, C, Ott, E. and Yorke, J. A., Chaotic attractors in crisis, Phys. Rev. Lett., 1982, 
48, 1507. 



29 



Guckenheimer, J. and Holmes, P., Nonlinear Oscillations, Dynamical Systems and Bifurca- 
tions of Vector Fields, 1986 (Springer- Verlag, New- York). 

Hughes D. W. and Proctor, M.R.E., A low-order model of the shear instability of convection: 
chaos and the effect of noise, Nonlinearity, 1990, 3, 127-153. 

Krishnamurti, R. and Howard, L. N., Large-scale flow generation in turbulent convec- 
tion.Proc. Natl. Sci. USA, 1981, 78,1981-1985. 

Labbe, R., Pinton, J. F. and Fauve, S., Study of the von Karman flow between coaxial 
corotating disks, Phys. Fluids, 1996, 8914-922. 

Larmor, J., How could a rotating body such as the sun become a magnet?. Rep. 87**^ Meeting 
Brit. Assoc. Adv. Sci., Bornemouth, Sept. 9-13, 1919, pp. 159-160, 1919 (London: John 
Murray) . 

Li, J., Sato, T. and Kageyama, A., Repeated and sudden reversals of the dipole fleld gener- 
ated by a spherical dynamo action. Science, 2002, 295, 1887-1890. 

Liu, B. and Zhang, J., Self-induced cyclic reorganization of free bodies through thermal 
convection, Phys. Rev. Lett, 2008, 100, 244501. 

Lorenz, E., Deterministic non periodic flow. Journal of the Atmospheric Sciences, 1963, 20, 
130-141. 

Mishra, P. K., De, A. K., Verma, M. K. and Eswaran, V., Dynamics of reorientations and 
reversals of large-scale flow in Rayleigh-Benard convection, J. Fluid Mech., 2010, (in press). 

Moffatt, H. K., Magnetic Field Generation in Electrically Conducting Fluids, 1978 (Cam- 
bridge University Press). 

Monchaux, R., Berhanu, M., Bourgoin, M., Moulin, M., Odier, Ph., Pinton, J.-F., Volk, R., 
Fauve, S., Mordant, N., Petrelis, F., Chiffaudel, A., Daviaud, F., DubruUe, B., Gasquet, C. 
and Marie, L., Generation of magnetic fleld by dynamo action in a turbulent flow of liquid 

sodium, Phys. Rev. Lett, 2007, 98, 044502. 

Mutabazi, I., Hegseth, J. J. and Andereck, C. D., Pattern formation in the flow between two 
horizontal coaxial cylinders with a partially flUed gap, Phys. Rev. A, 1988, 38, 4752-4760. 

Nishikawa, N. and Kusano, K., Simulation study of the symmetry- breaking instability and 
the dipole fleld reversal in a rotating spherical shell dynamo. Physics of Plasma, 2008, 15, 
082903. 



30 



Nozieres, P., Reversals of the Earth's magnetic field: an attempt at a relaxation model, 
Physics of the Earth and Planetary Interior, 1978, 17, 55-74. 

Paul, S, Kumar, K., Verma, M. K., Carati, D., De, A. K. and Eswaran, V., Chaotic travelling 
rolls in Rayleigh-Benard convection, Pramana J. Phys., 2010, 74, 75-82. 

Petrehs F. and Fauve S., Chaotic dynamics of the magnetic field generated by dynamo action 
in a turbulent fiow, J. Phys.: Condens. Matter, 2008, 20 494203. 

Petrelis, F. , Dormy, E., Valet, J.-P. & Fauve, S., Simple mechanism for the reversals of 
Earth magnetic field, 2009, Phys. Rev. Lett. 102, 144503. 

Petrelis F. and Fauve S., Mechanisms for magnetic field reversals. Philosophical Transactions 
of the Royal society A, 2010, 368, 1595-1605. 

Rabaud, M, Michalland, S. and Couder, Y., Dynamical regimes of directional viscous fin- 
gering: Spatiotemporal chaos and wave propagation, Phys. Rev. Lett., 1990, 64, 184-187. 

Ravelet, F., Marie, L., Chiffaudel, A. and Daviaud, F., Multistability and Memory Effect 
in a Highly Turbulent Flow: Experimental Evidence for a Global Bifurcation, Phys. Rev. 
Lett, 2004, 93, 164501. 

Ravelet, F., Berhanu, M., Monchaux, R., Aumaitre S., Chiffaudel, A., Daviaud, F, Dubrulle, 
B., Bourgoin, M., Odier, Ph., Plihon, N., Pinton, J.-F, Volk, R., Fauve, S., Mordant, N. and 
Petrelis F., Chaotic dynamos generated by a turbulent fiow of liquid sodium, Phys. Rev. 
Lett, 2008, 101, 074502. 

Rikitake, T., Oscillations of a system of disc dynamos, Proc. Camh. Phil. Soc, 1958, 54, 
89-105. 

Sarson, G.R. and Jones, C.A., A convection driven geodynamo reversal model. Physics of 
the Earth and Planetary Interior, 1999, 111, 3-20. 

Siggers, J. H., Dynamics of target patterns in low-Prandtl-number convection, J. Fluid 
Meek., 2003, 475, 357-375. 

Simon, A. J., Bechhoefer, J. and Libchaber, A., Solitary Modes and the Eckhaus Instability 
in Directional Sohdification, Phys. Rev. Lett, 1988, 61, 2574-2577. 

Sommeria J., Experimental study of the two-dimensional inverse energy cascade in a squarre 
box, J. Fluid Mech., 1986, 170, 139-168. 

Stefani F. and Gerbeth G., Asymmetric Polarity Reversals, Bimodal Field Distribution, and 



31 



Coherence Resonance in a Spherically Symmetric Mean-Field Dynamo Model, Phys. Rev. 
Lett, 2005, 94, 184506. 

Stieglitz, R. and Miiller, U., Experimental demonstration of a homogeneous two-scale dy- 
namo, Phys. Fluids, 200, 13, 561-564. 

Sugiyama, K., Ni, R., Stevens, R. J. A. M., Chan, T. S., Zhou, S.-Q., Xi, H.-D., Sun, C, 
Grossmann, S., Xia, K. Q. and Lohse, D., Flow Reversals in Thermally Driven Turbulence, 
Phys. Rev. Lett, 2010, 105, 034503. 

Tabeling P., Perrin B., Fauve S., Instability of a Linear Array of Forced Vortices, Europhys. 
Lett, 1987, 3 (4), 459-465. 

Tritton, D. J. Physical Fluid Dynamics 1977 (van Nostrand Reinhold Company, New York 
1977). 

Tuckerman, L. S. and Barkley, D., Global bifurcation to traveling waves in axisymmetric 
convection, Phys. Rev. Lett, 1988, 61, 408-411. 

Valet, J.-P., Meynadier, L. & Guyodo, Y. 2005 Geomagnetic field strength and reversal rate 
over the past 2 MiUion years. Nature 435, 802-805. 

VaUis, G. K., El Nino: A Chaotic Dynamical System? Science, 1986, 232, 243-245. 

Wicht, J. and Olson, P., A detailed study of the polarity reversal mechanism in a numerical 
dynamo model. Geochemistry, Geophysics, Geosystems, 2004, 5, Q03H10. 



32 



